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SUMMARY 


An  aircraft  flying  over  water  produces  a  significant  underwater 
sound  field  in  the  water  directly  beneath  it.  Acoustic  ray  theory  can  be 
used  to  obtain  a  useful  approximate  description  of  the  underwater  sound 
field  in  terms  of  sound  pressure  levels.  Details  of  the  theory  are 
presented.  Computer  programs  which  implement  the  theory  are 
documented.  The  results  include  time  histories  of  acoustic  signal  angles 
of  arrival  at  the  underwater  receiver;  time  histories  of  sound  pressure 
levels  received  at  various  depths;  and  equipressure  contour  descriptions 
of  the  underwater  sound  field. 
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Section  1 


INTRODUCTION 

An  aircraft  flying  over  water  produces  a  significant  underwater 
sound  field  in  the  water  directly  beneath  it. 

Given  a  measurement  of  the  aircraft  noise  source  level  in  air, 
acoustic  ray  theory  can  be  used  to  obtain  a  useful  approximate  description 
of  the  underwater  sound  field  in  terms  of  sound  pressure  levels.  The 
purpose  of  this  note  is  to  document  a  set  of  computer  programs  which 
have  been  developed  to  characterize  underwater  sound  fields  produced 
by  airborne  point  sources  under  specific  circumstances. 

Section  2  of  this  note  summarizes  the  existing  acoustic  ray  theory 
of  sound  transmission  from  air  into  a  body  of  water.  In  Section  3,  under¬ 
water  sound  fields  due  to  airborne  point  sources  are  characterized  in 
general.  Section  4  describes  the  set  of  computer  programs  which  have 
been  developed.  Finally,  Section  5  provides  comparisons  between  the 
predictions  of  ray  theoi’y  and  experimental  measurements,  which  provide 
validation  of  the  theory. 
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Section  2 

SOUND  TRANSMISSION  FROM  AIR  INTO  BODY  OF  WATER 
2. 1  EXISTING  ACOUSTIC  RAY  THEORY 

The  existing  acoustic  ray  theory  of  sound  transmission  from 
an  airborne  point  source  across  a  flat  air-water  interface  into  a  body 

of  water  has  been  well  established  for  about  15  years,  due  to  work  by 

12 
Hudimac  and  Horton. 

Consider  the  general  acoustic  ray  geometry  shown  in  Figure 
2-1.  There  is  a  stationary  point  sound  source  in  the  air  at  a  height 
above  the  flat  surface  of  the  water.  In  a  given  frequency  band  the 
source  sound  pressure  level,  p^,  is  known  at  some  reference 
distance,  rQ. 

The  sound  pressure  level  at  a  point,  R,  in  the  water  is  to  be 
determined.  In  the  acoustic  ray  theory  description,  sound  energy  is 
characterized  as  propagating  along  a  set  of  sound  rays  which  are 
emitted  radially  from  the  point  source,  and  which  are  redirected  upon 
crossing  the  air-water  interface.  The  basic  assumption  is  that  sound 
energy  does  not  cross  rays. 

Sound  pressure  level  pR  at  underwater  point  R  may  be 
determined  in  three  steps.  First,  the  reduction  in  sound  pressure 
due  to  spherical  divergence  as  the  sound  energy  propagates  along  a 
given  ray  from  the  source  to  the  water  is  determined.  Second,  the 
change  in  sound  pressure  as  the  sound  energy  is  transmitted  across 
the  air-water  interface  is  determined.  Third,  the  reduction  of  sound 
pressure  due  to  geometrical  spreading  as  the  sound  propagates  along 
a  given  ray  from  the  interface  to  the  receiver  at  point  R  is  determined. 
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The  first  reduction  in  sound  intensity  which  occurs  as  the 
acoustic  energy  propagates  through  the  air  is  given  by 

I1  =  ^rj0  /ri  )  I0> 

where  Iq  is  the  source  intensity  at  reference  distance  Tq,  and  1^  is 
the  sound  intensity  on  the  air  side  of  the  air-water  interface,  at  a 
distance  r^  from  the  source,  as  shown  in  Figure  2-1.  The  corre¬ 
sponding  sound  pressure  amplitude  is 

Pf  =  (r0/ri)  Pq. 

The  second  change  in  sound  intensity  as  the  acoustic  energy  is 

transmitted  across  the  air-water  interface  is  a  strong  function  of  the 

acoustic  impedances  of  the  air  and  water,  and  is  also  to  some  extent 

a  function  of  the  angle  of  incidence  of  the  acoustic  energy  upon  the 

surface  of  the  water.  Two  factors  are  of  importance.  These  are  the 

fraction  of  the  incident  acoustic  power  which  is  transmitted  across  the 

interface,  and  the  change  in  the  cross-sectional  area  of  the  ray  bundle 
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as  it  crosses  the  interface.  As  derived  in  detail  by  Brekhovskikh 
for  the  case  of  plane  waves,  the  relationship  between  incident  sound 
intensity,  1^,  and  transmitted  sound  intensity,  Ig,  is  given  by 

T  _  4 P1C1 P2C2 

o  a1> 

<»2C2  +  BP1C1) 

and  the  corresponding  incident  and  transmitted  sound  pressure 
amplitudes  are  related  by 


2p2C2 


2  (p2c2  +  Bp1c1)  ^1’ 


where  o^c^  and  p2c2  are  the  acoustic  impedances  of  the  air  and  water, 
respectively,  and  where  b  =  cos  0?/cos  0^. 
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Angle  ©2  is  related  to  angle  0^  through  Snell's  law: 
sin  @2  -  <c2/c1)  sin  0^ 

Thus,  numerical  factor  B  is  a  function  only  of  0^  and  the  sound  speeds 

c.,  and  C2  in  air  and  water,  respectively.  Snell's  law  is  a  consequence 
of  the  basic  physical  boundary  condition  requirements  that  the  sound 
pressure  amplitude  and  the  normal  component  of  particle  velocity  are 
continuous  across  the  air-water  interface. 

The  third  reduction  in  sound  intensity  as  the  acoustic  energy 
propagates  through  the  water  is  given  by 

(r.cosOJ 

t  = - _ - _ _  X 

(h1+c2h2/B3c1)  (h1+c2h2/Bc1)  2 

where  1^  is  the  sound  intensity  at  the  underwater  receiver  point  R, 
and  where  h1  and  h2  are  the  vertical  distances  shown  in  Figure  2-1. 
The  corresponding  sound  pressure  amplitude  at  point  R  is  given  by 

P2c2ricos01 

PR  (h1+c2h2/B3c1)l/2  (h1+c2h2/Bc1773  ?2* 

To  obtain  expressions  for  the  received  sound  intensities  and 
pressure  amplitudes  as  functions  of  the  corresponding  source  values 
and  the  appropriate  physical  parameters,  the  sets  of  three  respective 
intensity  and  pressure  expressions  indicated  above  are  combined.  The 
overall  results  are: 


*r 


pfflTQCoaeJ 


(h^Cghg/B  c^  (hj+Cghg/Bc^  (l+BpjCj/pgC^ 
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2rQcos01  pQ 

R  (hi+c2h2/B3ci)1/2  (hj-fC^/BCj)1/2  (1+BpjCj/ P2c2) 

2. 2  PHYSICAL  DESCRIPTION 

To  aid  in  the  visualization  of  underwater  sound  fields  due  to 
airborne  point  sources,  a  physical  interpretation  of  the  abstract 
mathematical  expressions  summarized  in  the  previous  section  is  needed. 

In  the  air,  there  is  spherical  divergence  of  sound  from  the 
point  source.  Sound  intensity  diminishes  as  the  inverse  square  of 
distance  from  the  source.  Sound  pressure  amplitude  diminishes  as 
the  inverse  of  distance  from  the  source.  Thus,  both  intensity  and 
pressure  diminish  by  6  dB  per  distance  doubled,  and  by  20  dB  per 
decade  increase  in  distance  from  the  source. 

At  the  air-water  interface,  any  acoustic  energy  propagating 
along  a  ray  which  subtends  an  angle  greater  than  about  13  degrees  from 
the  vertical  is  totally  reflected  from  the  surface  of  the  water.  As  indicated 
later  in  Figure  3-5,  only  acoustic  energy  within  a  cone  having  a 
half-angle  equal  to  the  critical  angle  of  about  13  degrees  is  transmitted 
into  the  water.  Even  in  this  region  where  the  critical  cone  intersects 
the  interface,  only  a  small  fraction  of  the  incident  acoustic  intensity  is 
transmitted.  The  transmission  loss  in  acoustic  intensity  is  about  29  dB. 

However,  it  is  important  to  note  that  there  is  not  a  correspond¬ 
ing  diminution  of  sound  pressure  amplitude  as  a  result  of  crossing  the 
air-water  interface.  The  acoustic  impedance  of  water  is  much  greater 
than  that  of  air.  Thus,  an  in-phase  reflection  of  the  incident  sound 
pressure  wave  occurs.  The  total  sound  pressure  amplitude  at  any 
instant  at  a  point  in  the  air  just  above  the  air-water  interface  is  almost 
double  the  instantaneous  amplitude  of  the  incident  sound  pressure  wave. 
Since  there  is  continuity  of  sound  pressure  across  the  air-water 
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interface,  this  means  there  is  sound  pressure  doubling  in  the  water 
just  below  the  air-water  interface  as  well.  Rather  than  being  decreased 
by  29  dB  upon  crossing  the  air-water  interface,  as  is  the  sound  intensity, 
the  sound  pressure  amplitude  is  actually  increased  by  6  dB  as  the 
acoustic  energy  crosses  the  interface.  Thus,  the  sound  pressure  levels 
in  the  water  are  about  35  dB  higher  than  the  corresponding  sound 
intensity  levels  in  the  water. 

The  above  situation  is  typical  for  sound  transmission  from  a 
low  impedance  medium  (air)  to  a  high-impedance  medium  (water).  The 
acoustic  intensity  diminishes  significantly,  while  the  sound  pressure 
amplitude  is  essentially  unchanged.  This  is  not  the  case  for  sound 
transmission  from  water  to  air.  The  acoustic  energy  reflected  from  a 
water-air  interface  is  out-of-phase  with  the  incident  energy.  Sound 
pressure  cancellation  occurs  in  the  water,  rather  than  pressure  doubling. 
Depending  upon  the  angle  of  incidence  at  the  water-air  interface,  the 
sound  intensity  diminishes  by  29  dB  or  more  mid  the  sound  pressure 
amplitude  diminishes  by  66  dB  or  more,  as  acoustic  energy  is  trans¬ 
mitted  from  water  into  air. 

Returning  now  to  die  case  of  sound  transmission  from  air  to 
water,  after  sound  crosses  the  air-water  interface,  its  intensity  and 
pressure  amplitude  diminish  further  as  it  geometrically  spreads  during 
propagation  through  the  vater.  Acoustic  energy  propagating  vertically 
downward  from  the  airborne  source  spreads  in  the  water  as  if  it  were 
radiated  from  a  point  on  the  vertical  axis,  at  a  height  on  the  vertical 
axis  given  by 


Acoustic  energy  entering  the  water  along  the  critical  ray,  however, 
spreads  in  the  vertical  plane  as  if  it  originated  from  a  source  which  was 
at  the  air-water  interface,  at  the  point  of  intersection  of  the  critical 


ray.  Thus,  most  of  the  acoustic  power  transmitted  across  the  air-water 
interface  propagates  nearly  vertically  downward.  The  amount  of  acoustic 
energy  redirected  in  near-horizontal  directions  is  relatively  small. 

2.3  SHORTCOMINGS  OF  EXISTING  THEORY 

The  principal  shortcoming  of  the  existing  acoustic  ray  theory 
of  underwater  sound  due  to  an  airborne  point  source,  as  developed  by 
Hudimac  and  by  Horton,  is  a  parametric  one.  An  expression  for  the 
received  sound  intensity  or  pressure  amplitude  is  available  only  in  terms 
of  the  source  values  and  the  following  physical  parameters: 
source  height,  hp  receiver  depth,  hg;  sound  speeds,  Cj,  and  c^; 
densities,  p^  and  p^;  and  the  angle  of  incidence,  0p  of  the  airborne 
sound  upon  the  air-water  interface. 

The  difficulty  is  that  the  horizontal  displacement  of  the 
receiver  from  a  vertical  line  through  the  source,  is  not  one  of  the  para¬ 
meters.  In  the  existing  form  of  the  theory,  one  cannot  specify  the 
source-interface-receiver  geometry  by  means  of  distances  hp  hg  and 
1 2  and  then  solve  directly  for  the  sound  intensity  or  pressure  amplitude 
at  the  receiver. 

Instead,  for  a  given  receiver  location,  it  is  necessary  to  use 

parameters  hp  hg  and  along  with  the  sound  speeds  c^  and  Cg  in 

Snell's  law  to  determine  the  angle  of  incidence,  0. .  This  requires  the 

2  1 

solution  of  a  quartic  equation  in  sin  0p  so  the  overall  task  is  not 
easily  accomplished  in  closed  form  by  hand  for  many  points. 

A  second  shortcoming  of  the  existing  theory  is  the  presentation 
of  the  results  in  terms  of  acoustic  intensity  only.  Sound  pressure  levels 
are  preferable.  In  experimental  acoustics,  sound  pressure  rather  than 
sound  intensity  is  the  practical  measurable  quantity.  In  air  to  water 
sound  transmission,  the  sound  pressure  and  intensity  levels  are  the 
same  in  air,  but  in  the  water  the  pressure  levels  are  about  35  dB  higher 
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than  the  intensity  levels,  due  to  the  change  in  acoustic  impedances. 

This  35  dB  difference  is  significant  when  underwater  sound  fields  due 
to  airborne  noise  sources  are  compared  with  underwater  sound  pressure 
levels  due  to  various  underwater  sources. 
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Section  3 


CHARACTERIZATION  OF  UNDERWATER  SOUND  HELD 
3. 1  EQUIPRESSURE  CONTOURS 


As  indicated  at  the  end  of  Section  2. 1,  the  received  underwater 
sound  pressure  amplitude  due  to  an  airborne  source  is,  according  to  acoustic 
ray  theory, 

2rp  cos  0i  Po 

PR  (h,  +  c2h2/B3c1)1/2  (h1+  c2h2/Bc1)1/2  (i  +  Bp  -^Ci/ P2c2) 


where 


and 


cos  @2 

B  '= - —  and  where 

cos 

Pq  =  source  sound  pressure  amplitude  at  r0; 
i'q  =  reference  distance  for  source  amplitude; 
0^  =  angle  of  ray  from  source  to  interface; 

@2  =  angle  of  ray  from  interface  to  receiver; 
hi  =  source  height  above  interface; 
h2  =  receiver  depth  below  interface; 
c^  =  sound  speed  in  air; 

C2  =  sound  speed  in  water; 

p  ^  =  density  of  air; 

p2  =  density  of  water; 

p^  =  received  sound  pressure  amplitude. 
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For  a  given  sound  source  and  air-water  medium,  parameters  Pq,  Tq, 
cl>  c2>  Pi  and  °2  have  constant  values.  The  three  basic  distances  required 
to  locate  an  underwater  receiver  point  relative  to  the  sound  source  are  height 
hj,  depth  h 2,  and  horizontal  offset  as  shown  in  Figure  3-1.  However, 
while  the  above  expression  for  p^  involves  h^  and  I12  directly,  it  is  not  an 
explicit  function  of  Distance  enters  only  implicitly,  through  angles 
and  6 2" 

To  evaluate  the  above  expression  for  p^  for  a  given  receiver  loca¬ 
tion,  that  is,  for  known  values  of  parameters  h^,  h2  and  it  is  necessary 
to  solve  for  0^  and  @2*  This  may  be  done  by  considering  the  geometry  of 
Figure  3-1,  and  by  using  Snell's  law  to  relate  @2  to  ®l  through  the  sound 
speeds,  C2  and  c^.  For  known  values  of  h^,  h2  and  the  values  of  0^, 

@2  and  are  unknown.  Three  equations  which  can  be  used  to  determine 
the  unknowns  are,  from  the  figure  and  from  Snell’s  law, 

j =  hj_  tan  6j; 

U2  -  j&l)  =  h2  tan  02; 

sin  02  =  ^/c^)  sin  el- 

The  problem  of  solving  the  above  three  equations  simultaneously  re- 

2 

duces  to  the  problem  of  solving  a  quartic  algebraic  equation  in  sin  0^: 

a^x4  +  agx3  +  a2x2  +  a^x  +  aQ  =  0 
where  x  s  sin2  0^  and  where 

a4  =  k4  [£24  +  {hj 2  -  h22)2  +  2/,22  (l^2  +  h22)]; 

a^  =  -2k3  [(k2  +  1)  +  k  h2 

+  (k2  +  2)  *  2hj2  +  (2k2  +  !)  *22h22 

-  (k2  +  1)  h  2h221; 
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where 


a2  -  [(k4  +  4k2  +  1)  ^24  +  hj4  +  k4h24 

+  2k2  (2k2  +  1)  l22h12  +  2k2  (k2  +  2)  j>22h22 
-  2k2h12h22]; 

a-^  =  -2[(k2  +  1)  x2^  +  i22h^2  +  k2  j^22h22];  and 

a  =  t  4 

ao  *2  ’ 

2  2 

k  =  (c^/c^)  in  the  above. 


When  a  computer  is  used  to  calculate  the  values  of  pR  at  a  number 
of  underwater  points,  the  above  quartic  algebraic  equation  is  easily  solved 
by  numerical  techniques.  Newton's  method  is  used.  A  first  approximation 
to  0p  chosen  from  the  range  0  g  01  s  13  degrees,  is  refined  iteratively  to 
the  desired  precision.  From  this  refined  value  of  0^,  the  values  of  02  and 
B  are  calculated.  Then  p  is  evaluated  directly. 


The  results  of  evaluating  pR  at  a  number  of  rectangular  grid  points 
serves  to  determine  the  underwater  sound  field,  as  shown  in  Figure  3-2. 
However,  a  more  convenient  representation  of  the  underwater  sound  field  is 
obtained  by  computing  and  plotting  equipressure  contours,  as  shown  in  Fig¬ 
ure  3-3  for  the  same  sound  pressure  data  presented  in  Figure  3-2.  The 
equipressure  contours  are  determined  by  means  oi  interpolation  between 
the  sound  pressure  levels  at  the  rectangular  grid  points. 

The  underwater  sound  field  shown  in  Figures  3-2  and  3-3  is  the 

O 

same  as  that  shown  by  Horton  .  The  source  height  is  300  feet.  The  source 
sound  pressure  level  is  p^  =  0  dB  at  reference  distance  r^  =  3  feet. 

Sound  pressure  levels  are  shown  in  Figures  3-2  and  3-3,  whereas  Horton 
shows  sound  intensity  levels,  which  are  about  35  dB  lower  because  of  the 
change  in  acoustic  impedance,  as  discussed  in  Section  2.  2. 
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3.2 


SOUND  PRESSURE  LEVELS  VERSUS  TIME 


The  underwater  sound  field  shown  in  Figures  3-z  and  3-3  has  rota¬ 
tional  symmetry  about  the  vertical  axis  through  the  omnidirectional  point 
source.  The  figures  show  vertical  slices  through  the  underwater  sound  field. 
The  vertical  axis  is  at  the  left,  so  that  horizontal  range  increases  to  the  right. 

In  the  majority  of  naval  applications  of  interest,  the  airborne  sound 
sources  are  airplanes  or  helicopters  which  move  horizontally  at  a  high  speed 
relative  to  the  underwater  receiver.  Time  histories  of  the  sound  pressure 
levels  at  the  underwater  receiver  are  useful  in  these  cases.  The  generation 
of  such  time  histories  may  be  visualized  by  moving  the  underwater  sound 
field  shown  in  Figure  3-3  past  an  underwater  receiver  which  is  located  at  a 
fixed  depth.  Alternatively,  the  receiver  may  be  visualized  as  moving  through 
the  fixed  sound  field  at  a  constant  depth,  at  the  speed  of  the  relative  hori¬ 
zontal  motion  between  the  source  and  receiver.  The  two  concepts  are  equiva¬ 
lent,  except  for  doppler  effects,  which  have  not  been  considered  in  the  present 
elementary  analysis. 

Figure  3-4  shows  some  typical  time  histories  of  received  underwater 
sound  pressure  levels  for  a  case  in  which  an  aircraft  flies  a  straight  hori¬ 
zontal  path  which  passes  directly  over  the  underwater  receiver  at  time  t  =  0. 
For  this  case,  the  aircraft  altitude  is  500  feet,  rather  than  300  feet  as  for  the 
case  shown  in  Figures  3-2  and  3-3.  Sound  pressure  level  time  histories  are 
shown  for  receiver  depths  of  100,  200,  400  and  800  feet.  At  the  greater 
depths,  the  maximum  received  sound  pressure  levels  diminish,  but  the  time 
intervals  over  which  the  time-varying  levels  exceed  some  arbitrary  threshold 
level  are  increased  in  duration. 

The  case  in  which  the  aircraft  flies  directly  over  the  underwater  re¬ 
ceiver  is  an  extreme  example.  More  often,  there  will  be  a  horizontal  offset 
between  the  path  of  the  aircraft  and  the  vertical  axis  through  the  array.  Time 
histories  of  the  received  sound  pressure  levels  for  the  case  of  a  nonzero  offset 
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of  the  source  are  just  as  easily  calculated  as  for  the  case  of  zero  offset. 
However,  to  visualize  the  generation  of  the  time  histories  of  received  sound 
pressure  levels  in  the  case  of  nonzero  offset,  it  is  necessary  to  consider  the 
three-dimensional  underwater  sound  field  rather  than  a  two-dimensional 
vertical  radial  slice  through  it. 

3.  3  INCLUSION  OF  RECEIVER  NOISE  EFFECTS 

In  most  of  the  applications  of  interest,  the  relationship  between  the 
transmitted  underwater  noise  signal  of  the  aircraft  and  the  self-noise  of  the 
underwater  receiver  is  crucial.  Generally,  the  total  signal-plus-noise  output 
from  the  receiver  will  be  frequency  band-limited,  so  as  to  concentrate  on 
spectral  regions  where  the  aircraft  radiated  noise  signal  is  strong.  Then, 
two  parameters  are  significant.  These  are  the  peak  signal-plus-noise  to 
noise  ratio,  and  the  duration  of  the  time  interval  over  which  the  signal-plus- 
noise  level  exceeds  the  noise  level  by  some  preselected  amount,  such  as  3  dB. 

The  self-noise  of  the  receiver  may  be  flow  noise,  in  the  case  of  hull- 
mounted  or  towed  underwater  receivers,  or  it  may  be  the  ambient  ocean  noise 
in  the  frequency  band  of  interest.  In  either  case,  it  may  be  safely  assumed 
that  the  noise  adds  to  the  aircraft  signal  incoherently.  This  means  that  the 
noise  and  signal  powers  are  added  before  conversion  to  decibels.  The  re¬ 
sults  are  signal-plus-noise  level  time  histories  which  are  similar  to  those 
shown  in  Figure  3-4,  but  which  are  limited  at  their  lower  bounds  by  the  con¬ 
stant  noise  level  which  has  been  chosen. 

3.4  GEOMETRY  OF  ACOUSTIC  WAVEFRONTS 

The  geometry  of  the  acoustic  wavefronts  involved  in  the  generation 
of  an  underwater  sound  field  by  an  airborne  point  source  is  shown  in  Fig¬ 
ure  3-5.  This  diagram  was  constructed  using  only  acoustic  ray  theory. 

In  the  diagram,  acoustic  rays  are  shown  as  straight-line  segments, 
uniformly  distributed  at  5-degree  increments  in  the  water.  Snell's  law  was 
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used  to  determine  the  bending  of  the  rays  at  the  air-water  interface.  Those 
rays  along  which  sound  is  transmitted  from  the  air  into  the  water  are  en¬ 
closed  in  a  vertical  cone,  the  half-angle  of  which  is  approximately  13  de¬ 
grees.  The  region  of  sound  transmission  is  indicated  in  the  diagram  by  an 
encircled  "T.  "  As  discussed  in  an  earlier  section,  most  of  the  incident 
sound  energy  inside  the  cone  of  transmission  is  reflected  back  upward  at  the 
air-water  interface,  and  yet  the  sound  pressure  amplitude  (not  energy  or  in¬ 
tensity)  transmitted  across  the  interface  is  appreciable,  comparable  in  mag¬ 
nitude  to  the  incident  sound  pressure  amplitude. 

Outside  the  cone  of  transmission,  the  acoustic  energy  incident  upon 
the  air-water  interface  is  totally  reflected.  Only  a  slight  penetration^  of 
sound  into  the  water  occurs.  This  penetrating  sound  dies  out  exponentially 
with  depth  and  with  frequency;  this  is  not  a  sound  propagation  effect-.  The 
region  of  sound  penetration  is  indicated  in  the  diagram  by  an  encircled  "P.  " 

Acoustic  wavefronts  are  shown  in  Figure  3-5  as  curved  segments 
which  are  perpendicular  to  the  rays.  The  wavefronts  are  shown  at  equally- 
spaced  successive  instants  of  time.  One  unit  of  time  is  the  interval  between 
the  time  at  which  the  spherically  spreading  airborne  acoustic  wavefront  first 
reaches  the  water  directly  below  the  source,  and  the  time  at  which  the  wave- 
front  impinges  upon  the  water  at  the  critical  angle  of  approximately  13  de¬ 
grees.  This  duration  is  directly  proportional  to  the  height  of  the  source 
above  the  water.  For  any  source  height,  as  shown  in  the  diagram  it  requires 
about  35.  5  such  intervals  of  time  for  sound  to  propagate  from  the  source  to 
the  water  directly  below.  Thus,  one  unit  of  time,  At,  in  Figure  3-5  is  given 
approximately  by 


where  h^  is  the  source  height  and  c^  is  the  speed  of  sound  in  air. 
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3.5 


ANGIE  OF  ARRIVAL  OF  SIGNAL  AT  RECEIVER 


As  in  Section  3. 1,  let  the  position  of  the  underwater  receiver  rela¬ 
tive  to  the  airborne  source  be  described  in  terms  of  source  height  hp  re¬ 
ceiver  depth  hg,  and  horizontal  offset  Then,  as  described  earlier,  a 
quartic  algebraic  equation  in  sin^  0^  can  be  solved  to  determine  0p  the 
pola>-  angle  of  arrival  of  the  signal  at  the  underwater  receiver.  This  is  an 
angle  measured  relative  to  a  vertical  axis  through  the  receiver. 

Usually,  the  underwater  receiver  has  a  definite  horizontal  axis  of 
symmetry.  This  axis  will  have  a  certain  azimuthal  heading,  relative  to  the 
heading  along  which  the  airborne  sound  source  moves.  The  angle  of  arrival 
of  the  acoustic  signal  relative  to  the  horizontal  axis  of  the  receiver  can  then 
be  determined  from  the  polar  angle,  0p  and  the  relative  heading,  8,  of  the 
receiver,  by  consideration  of  the  solid  geometry  which  is  involved. 

Specifically,  suppose  that  the  underwater  receiver  is  a  horizontal 
line  array  of  hydrophones.  The  axis  of  symmetry  of  this  receiver  is  the 
horizontal  line  along  which  the  hydrophones  lie.  With  such  an  array,  it  is 
possible  to  measure  only  the  angle  of  arrival  of  an  acoustic  signal  relative  to 
the  axis  of  the  array.  It  is  not  possible  to  measure  the  vertical  angle  of  ar¬ 
rival  of  the  signal  at  the  array. 

As  a  convention,  let  the  heading  of  the  underwater  receiver  relative 
to  the  flight  path  of  the  aircraft  be  zero  when  the  direction  of  the  axis  of  the 
receiver  is  parallel  and  opposite  to  the  direction  of  flight  of  the  aircraft. 

With  this  convention,  when  the  flight  path  passes  directly  over  the  receiver, 
with  the  receiver  heading  zero  degrees,  the  angle  of  arrival  of  the  acoustic 
signal  relative  to  the  axis  of  the  array,  0^,  varies  from  an  initial  value  of 
zero  to  a  final  value  of  180  degrees.  Similarly,  when  the  receiver  heading 
is  90  degrees,  measured  either  clockwise  or  counterclockwise  from  zero, 

0^  -  90  degrees  throughout  the  run,  when  the  flight  path  passes  directly 
over  the  receiver. 
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Figure  3-6  shows  time  histories  of  0 ^  for  a  flight  passing  directly 
over  the  receiver.  The  source  height  is  500  feet,  the  same  as  in  Figure  3-4 
where  sound  pressure  level  time  histories  were  shown.  The  receiver  depth 
in  Figure  3-6  is  400  feet.  Results  shown  in  the  figure  are  consistent  with 
the  preceding  discussion. 

Figure  3-7  shows  similar  angles  of  arrival  for  a  source  height  of 
4000  feet.  As  would  be  expected  on  the  basis  of  the  acoustic  ray  geometry 
shown  in  Figure  3-5,  when  the  source  is  at  a  greater  altitude  the  entire 
underwater  sound  field  pattern  scales  up  proportionally  larger.  At  a  depth 
of  400  feet,  the  angles  of  arrival  vary  more  slowly  with  time,  for  the  same 
aircraft  speed. 

A  more  interesting  case  is  shown  in  Figure  3-8.  Here,  the  condi¬ 
tions  are  the  same  as  for  Figure  3-6,  except  that  the  flight  path  does  not 
pass  directly  over  the  underwater  receiver.  The  horizontal  offset  of  the 
flight  path  from  the  vertical  axis  through  the  receiver  is  2000  feet  at  the 
closest  point  of  approach.  As  a  result,  the  time  histories  of  angles  of  ar¬ 
rival  of  the  signal  in  Figure  3-8  are  somewhat  altered  from  those  in  Fig¬ 
ure  3-6.  The  angle  of  arrival  is  not  constant  for  any  receiver  heading,  when 
there  is  a  nonzero  offset  of  the  flight  path.  This  case  is  typical.  The  case 
of  a  flight  directly  over  the  underwater  receiver  is  an  extreme  example. 

3.6  LIMITING  CASES  FOR  UNDERWATER  SOUND  HELD 


As  summarized  in  detail  at  the  start  of  Section  3. 1,  the  acoustic 
ray  theory  expression  for  the  received  sound  pressure  amplitude  due  to  an 
airborne  point  source  is 


2rQ  cos 


PR  ~  9  i  /o  1  /o 

(hj  +  c2h2/B  c^  (hj  +  c^/BCj)  7  (1  +  B^c^p^) 


in  which  B  =  cos  02/cos  0^  and  where  the  other  parameters  are  defined  in 
Section  3.1.  This  expression  has  been  used  in  the  computer  programs 
described  in  this  note. 
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ANGLE  OF  ARRIVAL 
AT  RECEIVER  AXIS 
(deg) 


Histories  of  Angles  of  Arrival  of  Signal  at  Underwater  Receiver 


ECEIVER  AXIS 


Consider  now  some  basic  physical  approximations  in  the  preceding 
equation.  Within  the  cone  of  transmission,  the  angle  0^  takes  on  values 
0  £  0i  €  13  degrees.  Correspondingly,  cos  0^  takes  on  values 
1  -  cos  8^  =  0.97,  so  for  all  practical  purposes  cos  6-^  may  be  approxi¬ 
mated  as  unity,  and  B  =  cos  02  is  a  good  approximation.  Thus,  dimen¬ 
sionless  factor  B  takes  on  values  0  <  B  g  1. 


Also,  for  air  and  water  the  ratio  p^ / P2C2  has  a  value  of  about 
0.0003.  Thus,  the  approximation  (1  +  Bp  1C1/P2C2)  =  1  is  a  good  one. 

Using  the  latter  approximation,  and  with  the  range  of  values  of  B 
in  mind,  a  good  approximation  to  p^  is  the  expression 


PR  " 


2BVo 


(B3h1  +  c2h2/c1)1/2  (Bh^Cghg/c^172 


This  approximation  is  valid  for  all  permissible  values  of  p^,  r^,  h^, 
and  B. 


For  the  first  limiting  case,  consider  the  situation  where  the  re¬ 
ceiver  is  located  directly  under  the  source.  In  this  case,  0^  =  0,  @2  =  0, 
and  B  =  1.  The  sound  pressure  amplitude  received  in  the  water  directly  be¬ 
low  the  source  is 


2ropo 


2ropo 


pR  "  (hj  +  (hj  +  5h2)’ 


since  the  ratio  Cg/c^  is  approximately  5  for  water  to  air.  . 

If  the  receiver  depth  is  small  compared  with  the  source  height,  that 
is,  if  h^  ^  h^/50,  then  the  received  sound  pressure  amplitude  is  a  function 
of  the  source  height,  h^,  only: 

r0 

pR  =  2  yj-  Pq  when  h2  -  h^/50. 
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This  sound  pressure  amplitude  just  under  the  surface  of  the  water,  directly 
below  the  sound  source,  is  twice  the  amplitude  in  the  air  just  above  the  sur¬ 
face  of  the  water,  due  to  pressure  doubling  at  the  interface. 

If  the  receiver  depth  is  greater  than  double  the  source  height,  then 
the  received  sound  pressure  amplitude  is  a  function  of  the  receiver  depth,  h^, 
only: 

r 

pR  =  0. 4~P0  when  h2  2  2hr 

So  long  as  the  receiver  depth  is  greater  than  double  the  source  height,  the 
underwater  sound  field  pressure  levels  diminish  by  6  dB  per  depth  doubled, 
directly  below  the  sound  source. 

Thus,  directly  under  the  sound  source,  the  sound  pressure  levels 
just  under  the  surface  of  the  water  are  constant  (at  double  the  value  in  the 
air  just  above  the  surface  of  the  water)  down  to  a  depth  of  roughly  h2  =  h]/50. 
The  sound  pressure  levels  then  diminish  with  increasing  depth,  at  a  rate  which 
increases  with  depth  from  only  a  fraction  of  one  dB  per  distance  doubled  at 
the  depth  h2  =  h^/50,  to  a  rate  of  6  dB  per  distance  doubled  at  the  depth 
h£  =  2hj.  At  depths  greater  than  h2  =  2hj,  directly  under  the  sound 
source,  the  sound  pressure  levels  diminish  with  depth  at  a  constant  rate  of 
6  dB  per  distance  doubled. 

For  the  second  and  more  general  limiting  case,  consider  situations 
in  which  the  receiver  is  located  far  from  the  source,  in  units  of  source  height. 
Referring  to  the  earlier  Figure  3-1,  let  r2  be  the  slant  distance  from  the 
point  where  the  acoustic  ray  crosses  the  air-water  interface,  to  the  point 
where  the  receiver  is  located.  That  is,  let  r2  =  h2/cos  e2.  If  this  distance 
r2  is  greater  than  double  the  source  height,  then  the  expression  for  the  re¬ 
ceived  sound  pressure  amplitude,  pR,  can  be  simplified. 
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This  can  be  seen  by  consideration  of  the  case  in  which  Bh^  «  C2h2/c^ 
If  this  condition  is  satisfied,  then  the  condition  B^h^  <<  C2h2/c-^  is  also  satis¬ 
fied,  since  0  <  B  S  1,  and  the  received  sound  pressure  amplitude  becomes 


PR  " 


25  Vo 

C2h2/Cl 


=  0.4B 


2^0 

h2  °’ 


where  c2/c^  =  5  is  assumed. 

Since  B  =  cos  @2  =  h2//r2’  receivec*  sound  pressure  amplitude 
for  this  case  becomes 


r0h2 

pR  =  0.4-2/  p 


0* 


The  condition  under  which  the  above  expression  is  valid  is  Bh.  -c  c  h  /c 

1  Li  Li  JL 

or  approximately  r2  s  2hr  That  is,  when  the  slant  distance  r2  from  the 
surface  below  the  source  to  the  receiver  is  greater  than  double  the  source 
height,  hp  the  above  expression  is  valid.  The  above  expression  indicates 
that  at  distances  rg  £  2h^,  the  received  sound  pressure  level  at  a  given 
depth,  h2,  diminishes  at  the  rate  of  12  dB  per  distance  (r2)  doubled. 

Directly  below  the  sound  source,  r2  =  h2,  and  the  above  expres¬ 
sion  for  pR  at  r2  g  2h^  becomes 


This  agrees  with  the  result  obtained  previously,  for  the  case  directly  below 
the  source,  at  a  depth  h2  s  2hj. 

In  general,  it  can  be  shown  that  the  equipressure  contours  which  de¬ 
scribe  the  underwater  sound  field  due  to  an  airborne  point  source  take  the 
limiting  form  of  circles,  in  the  more  distant  regions  where  r2  i!  2h1.  The 
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circles  are  all  tangent  to  the  surface  of  the  water,  directly  below  the  sound 
source.  The  radii  of  the  circles  are  easily  determined  from  the  values 

PR  =  0.4^  p0 

of  the  sound  pressure  amplitudes  on  the  vertical  axis  passing  through  the 
source. 

Figure  3-9  is  an  example  of  an  actually  computed  underwater  sound 
field  in  the  limiting  case  in.  which  the  airborne  sound  source  is  close  to  the 
surface  of  the  water.  The  equipressure  contours  are  clearly  circles  which 
are  tangent  at  the  air-water  interface,  directly  below  the  source.  In  this 
case,  the  source  was  assumed  to  be  a  6- grain  dynamite  cap  which  was  ex¬ 
ploded  at  a  height  of  one  foot  above  the  water.  The  purpose  of  computing 
this  particular  underwater  sound  field  was  to  investigate  the  possibility  of 
using  near- surface  airborne  explosive  charges  as  underwater  transient  sound 
sources.  Such  a  technique  could  avoid  surface  reflections,  bubble  pulses  and 
fish  kills,  all  of  which  are  problems  with  underwater  explosions.  Also,  this 
technique  could  direct  the  sound  downward  and  would  limit  the  amplitude  of 
the  horizontally  propagating  portion  of  the  sound. 

The  expressions  obtained  in  the  above  limiting  cases  are  useful  in 
making  approximate  calculations  of  underwater  sound  fields  due  to  airborne 
sources.  However,  the  exact  expressions,  such  as  that  for  pR  given  at  the 
beginning  of  this  section,  have  always  been  used  in  the  sequence  of  computer 
programs  described  in  this  note. 


3-20 


PRECEDING  PAGE  BLANK-NOT  FILMED 


Section  4 

SEQUENCE  OF  COMPUTER  PROGRAMS 

A  sequence  of  seven  computer  programs  has  been  developed 
to  compute  and  machine  plot  various  functions  associated  with  under¬ 
water  sound  fields  generated  by  airborne  point  sources.  The  programs 
have  been  written  in  Super  Basic  for  use  on  a  time-sharing  computer 
terminal. 

Detailed  listings  of  the  seven  computer  programs  are  given 
in  Appendix  A.  The  purpose  of  the  present  section  is  to  describe  the 
programs  from  the  user's  point  of  view. 

Program  USF  (Underwater  Sound  Field)  is  the  basic  program 
in  the  sequence.  It  is  used  to  compute  received  sound  pressure  levels 
in  the  water  at  points  on  a  rectangular  grid,  as  shown  previously  in 
Figure  3-2.  Vertical  angles  of  arrival  of  the  signal  at  the  various 
points  on  the  grid  are  determined  as  a  part  of  the  computation.  Computed 
sets  of  values  are  outputted  in  data  files,  ready  for  input  into  successive 
programs. 

Program  AA  (Angle  of  Arrival)  is  used  to  compute  time 
histories  of  angles  of  arrival,  0A,  of  the  signal  relative  to  the  horizontal 
axis  of  the  receiver,  for  various  receiver  depths  and  headings.  The 
values  of  the  angles  are  outputted  to  a  succeeding  program,  APLOT,  for 
machine  plotting  of  the  results.  Typical  results  have  been  shown  r 
Figure  3-6,  3-7  and  3-8. 

Program  ADDAMB  (Add  Ambient)  makes  possible  the  inclusion 
of  ambient  noise  levels  in  the  underwater  sound  rield  data  computed  in 
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program  USF.  The  ambient  noise  level  which  is  selected  is  combined 
with  the  noise  levels  due  to  the  airborne  source  in  a  sound  power 
addition,  prior  to  decibel  conversion. 

Program  LPLOT  {Levels  Plot)  is  used  to  prepare  time 
histories  of  received  sound  pressure  levels  at  various  depths,  from 
program  USF,  modified  by  inclusion  of  the  ambient,  for  machine 
plotting.  Typical  results  have  been  shown  in  Figure  3-4. 

Finally,  program  CON  (Contours)  is  used  to  compute  under¬ 
water  sound  field  equipressure  contours  from  the  sound  field  data  from 
program  USF,  modified  by  the  inclusion  of  the  ambient.  An  interpolative 
process  is  used.  The  output  from  CON  is  prepared  for  machine  plotting 
by  program  CPLOT.  Typical  results  have  been  shown  in  Figure  3-3. 

The  three  plotting  programs,  APLOT,  CPLOT  and  LPLOT, 
are  used  to  properly  scale  the  data,  assign  the  plot  origin,  and  to  carry 
out  other  similar  functions  in  preparing  the  data  for  machine  plotting. 

The  data  to  be  plotted  is  outputted  from  data  files  on  paper  tape  at  the 
time- sharing  computer  terminal.  It  is  then  inputted  to  a  minicomputer 
system  which  is  operated  in-house  at  Hydrospace,  for  machine  plotting 
on  a  digital  incremental  plotter.  As  the  final  step  in  plot  preparation, 
machine  plots  of  data  are  pasted  onto  preprinted  sheets  having  the 
appropriate  ordinate  and  abscissa  scales. 

The  flow  diagram  for  the  sequence  of  seven  computer  programs 
is  shown  in  Figure  4-1.  The  input  parameters  are  different  for  the 
various  programs;  the  input  parameters  for  each  program  are  shown  in 
Tables  4-1  and  4-2.  The  input  and  output  data  files  for  each  of  the 
seven  programs  are  summarized  in  Table  4-3. 
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Figure  4-1.  Flow  Diagram  for  Sequence  of  Computer  Programs 
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Table  4-1.  Input  Parameters 


Parameter 


PI 

Rl 

HI 

D1 


U 


Y 


SI 

PO 


Units 


Description 


dB  //  ubar  Source  pressure  level  in  air  at 
distance  Rl 


ft  Reference  distance  at  which 

PI  is  measured 


ft  Height  of  source  above  water 

ft  Horizontal  offset  distance  from 

source  to  vertical  axis  through 
receiver  at  closest  point  of 
approach  of  aircraft  flight  path 

ft  Grid  increment  for  rectangular 

grid  of  points  at  which  under¬ 
water  sound  field  is  determined 


ft 


knots 

dB  //ubar 


Maximum  horizontal  distance  to 
which  underwater  sound  field  is 
determined  in  Program  USF 

Speed  of  aircraft 

Ambient  sound  pressure  level 


Y2 


C2 


Cl 


ft  Maximum  horizontal  distance  to 

which  equipressure  contours  are 
to  be  determined  in  Program  CON 

dB  Interval  between  equipressure 

contours  of  underwater  sound 
field  in  Program  CON 

dB  Interval  between  equipressure 

contours  of  airborne  sound  field 
in  Program  CPLOT 
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Table  4-2.  Parameters  Inputted  to  Each  Program 


Program 

Input  Parameters 

USF 

Pi,  Rl,  Hi,  Dl,  U,  Y 

AA 

(none) 

APLOT 

SI 

ADDAMB 

PO 

LPLOT 

SI 

CON 

Y2,  C2 

CP  LOT 

Cl 

Table  4-3. 

Input  and  Output  Data  Files 

Program 

Input  Data  Files 

Output  Data  Files 

USF 

(none) 

ADATA 

CDATA 

UYDATA 

AA 

UYDATA 

ADATA 

BDATA 

APLOT 

UYDATA 

BDATA 

APDATA 

ADDAMB 

UYDATA 

CDATA 

UYDATA 

FDATA 

LPLOT 

UYDATA 

FDATA 

LPDATA 

CON 

UYDATA 

FDATA 

UYDATA 

EDATA 

DDATA 

CP  LOT 

UYDATA 

DDATA 

EDATA 

CPDATA 
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It  is  anticipated  that  the  foregoing  rather  brief  description 
will  suffice  as  an  introduction  to  the  use  of  the  sequence  of  seven 
computer  programs  for  the  determination  of  underwater  sound  fields 
due  to  airborne  sources.  Each  of  the  computer  programs  is  extensively 
self-explanatory.  Upon  execution,  each  program  poses  to  the  user  a 
series  of  questions,  such  as,  ,fWhat  is  the  source  pressure  level,  PI, 
in  dB  //  ubar  at  Rl?".  These  questions  assure  the  efficient  and  error- 
free  inputting  of  the  input  variables.  The  best  way  of  gaining  further 
familiarity  with  the  sequence  of  computer  programs  is  to  put  them  to 
actual  use. 
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Section  5 


COMPARISON  BETWEEN  RAY  THEORY  PREDICTIONS 
AND  EXPERIMENTAL  MEASUREMENTS 

This  report  has  described  a  sequence  of  computer  programs 

which  implement  the  acoustic  ray  theory  predictions  developed  by  Hudimac'*- 
2 

and  by  Horton  for  underwater  sound  fields  due  to  airborne  sources. 

As  indicated  by  Hudimac\  the  ray  theory  in  this  case  will  not 
be  valid  for  large  horizontal  separations  {n^)  of  the  source  and  receiver, 
especially  when  the  receiver  is  near  the  surface.  That  is,  the  ray  theory 
predictions  may  be  inaccurate  when  the  receiver  is  located  within  an 
acoustic  wavelength  of  the  surface,  at  the  source  frequency  of  interest. 

In  these  cases,  acoustic  wave  theory  rather  than  ray  theory  should  be 
applied. 

However,  the  validity  of  the  acoustic  ray  theory  predictions  of 

the  underwater  sound  field  directly  beneath  the  aircraft  has  also  been 

4 

questioned  in  recent  years  by  Weinstein  and  Henney  .  In  this  work,  the 
zero- frequency  limit  of  the  wave  theory  expression  for  the  underwater  sound 
field  directly  beneath  the  aircraft  was  taken,  and  it  did  not  agree  with  the 
corresponding  ray  theory  expression. 

5 

Young  has  recently  made  some  experimental  acoustic  measure¬ 
ments  to  determine  which  is  valid,  Hudimac's  ray  theory  or  Weinstein's 
limiting  form  of  wave  theory,  for  the  prediction  of  sound  pressure  levels 
at  points  in  the  water  directly  below  an  airborne  point  source.  It  has  been 
concluded  that  ray  theory  provides  an  adequate  description  in  this  case, 
and  the  limiting  form  of  the  wave  theory  does  not.  It  appears  that  the  zero- 
frequency  limit  is  too  much  of  an  extreme  case  to  apply  to  noise  at  useful 
frequencies,  above  10  Hz. 
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A  further  comparison  between  ray  theory  predictions  and  experi¬ 
mental  measurements  has  been  made  recently  at  Hydrospace,  using  data 

g 

from  a  sonic  boom  simulation  experiment  in  which  the  penetration  of  air¬ 
borne  shock  wave  energy  into  a  body  of  water  was  the  principal  subject  of 
investigation.  (Here,  the  phenomenon  of  penetration  of  an  acoustic  dis¬ 
turbance  across  the  air-water  interface,  in  the  form  of  a  non- propagating 
evanescent  wave  accompanying  a  totally- reflected  airborne  wave,  is  dis¬ 
tinguished  from  the  phenomenon  of  transmission  of  sound  across  the  air- 
water  interface.)  In  this  work,  dynamite  caps  were  exploded  over  a  flooded 
quarry.  Transmission  effects  were  measured,  in  addition  to  penetration 
effects. 

Figure  5-1  shows  the  comparison  between  the  experimental  measure¬ 
ments  and  the  predictions  of  ray  theory,  as  determined  using  the  sequence 
of  computer  programs  described  in  this  note.  The  experimental  values 
are  enclosed  with  rectangles  and  are  plotted  at  the  points  of  measurement. 

The  theoretical  values  are  enclosed  with  ovals  and  are  shown  as  equipressure 
contours  of  the  underwater  sound  field.  The  source  height  was  30  ft. 

Agreement  between  the  experimental  measurements  and  the  pre¬ 
dictions  of  ray  theory  as  shown  in  Figure  5-1  is  quite  good,  within  +,  3  db 
in  most  cases  throughout  the  region  shown.  This  good  agreement  further 
validates  acoustic  ray  theory  for  the  prediction  of  underwater  sound  fields 
due  to  airborne  point  sources.  The  investigator  who  uses  the  sequence  of 
computer  programs  described  in  this  note  to  predict  underwater  sound 
fields,  given  the  altitude  and  source  level  of  the  airborne  noise  source, 
can  do  so  with  confidence  that  the  predictions  will  be  accurate  except  at 
points  which  are  very  distant  horizontally,  especially  close  to  the  air- 
water  interface. 


Figure  5-1.  Comparison  Between  Ray  Theory  Predictions  and  Experimental 
Measurements  of  Underwater  Sound  Field  Due  to  Airborne  Point 
Source  at  Height  of  30  feet 
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APPENDIX  A 


LISTINGS  0F  COMPUTER  PROGRAMS 
FOR  DETERMINATION  OF 
UNDERWATER  SOUND  FIELDS 
DUE  TO  AIRBORNE  POINT  SOURCES 
USING  ACOUSTIC  RAY  THEORY 


PROGRAMS 
LI  STED: 

/USF/ 

/AA/ 

/APLOT/ 

/ADD AM B/ 

/LPLOT/ 

/CON/ 

/CPL0T/ 


PROGRAMS  ARE  WRITTEN  IN 
"SUPER  BASIC" 
LANGUAGE  FOR  USE  ON 
"TYM  SHARE" 

TIME- SHARING  COMPUTER 


PROGRAM  /USF / 


1000  PRINT  "THIS  IS  PROGRAM  /USF/,  FOR  COMPUTING  THE  UNDERWATER  SOUND” 
1010  PRINT  "  FIELDS  ASSOCIATED  WITH  AIRBORNE  SOUND  SOURCES" 

1020  PRINT 

1030  PRINT  "DATE  AND  TIME:  ":  DATE 

1040  PRINT 

1050  PRINT 

1060  VAR  =  ZERO 

1070  G$  =  TIME  /60 

1080  DIM  L $( 0: 60) ,  FSC0:60>,  P$(  1 :  1 0* 0: 60> ,  T$(  1 :  1 0,  0:  60) 

1090  INTEGER  PI,  Rl,  HI,  Dl,  U,  Y,  Sl>  P0>  K9 

1100  VI  =  1085 

1110  V 2  =  4750 

1120  Rl  =  1 .29 

1130  R2  =  1025 

1  140  K1  =  V2  /  VI 

1150  Z  1  =  (  V 1* R 1 )  /  C  V2*R2) 

1160  K2  =  K 1 : 2 
1170  K 4  =  K2T2 

1180  PRINT  "WHAT  IS  SOURCE  PRESSURE  LEVEL,  PI,  IN  DB//UBAR  AT  Rl": 

119  0  INPUT  P! 

1200  PRINT 

1210  PRINT  "WHAT  IS  REFERENCE  DISTANCE,  Rl,  IN  FEET": 

1220  INPUT  Rl 
1230  PRINT 

1240  PRINT  "WHAT  IS  HEIGHT,  HI,  0F  SOURCE  ABOVE  WATER,  IN  FEET": 

1250  INPUT  HI 
1260  PRINT 

1270  PRINT  "WHAT  IS  OFFSET  DISTANCE,  Dl,  FROM  SOURCE  TO  VERTICAL  AXIS" 

1280  PRINT  "  THROUGH  RECEIVER  AT  CLOSEST  POINT  OF  APPROACH,  IN  FEET": 

1290  INPUT  Dl 

1300  Q 1  =  K2  +  1 

1310  Q2  =  K2  +  2 

1320  Q3  =  2  *  K2  +  1 

1330  Q4  =  K4  +  4  4.  K2  +  1 

1340  Q5  =  2  *  K2 

1350  Q6  =  2  *  03 

1360  Q7  =  Q5  *  92 

1370  PRINT 

1380  PRINT  "WHAT  IS  GRID  INCREMENT,  U,  IN  FEET": 

1390  INPUT  U 


A- 2 


1400  PRINT 

1410  PRINT  "WHAT  IS  MAXIMIM  DISTANCE*  Y*  IN  FEET”: 

1 4 20  INPUT  Y 
1430  PRINT 
1440  J9  =  Y  /  U 
1450  FOR  J  =  0  T0  J9 
1460  D2  =  J  *  U 

1470  IF  D1  =  0  AND  D2  =  0  THEN  FSCJ)  =  0 

1480  IF  D1  =  0  AND  D2  >  0  THEN  FSCJ)  =  90 

1 49 C  IF  D1  >  0  AND  D2  >  0  THEN  FSCJ)  *  C  180/PI)  *  ATANC  D2/D1  ) 

1500  LSCJ)  =  SQRTC  Dl»2  +  D2r2  ) 

1510  IF  LSCJ)  =  0  THEN  LSCJ)  =  1 E-  6 
1520  NEXT  J 

1530  P8  =  10  t  ( P 1/20) 

1540  F0R  I  =  1  T0  10 
1550  M2  -  I  *  U 
1550  PRINT 
«5 70  PRINT  *' 

DEPTH,  H2*  0F  RECEIVER*  IN  FEET*  = 

1530  PRINT  IN  IMAGE  ,,XXXX,,i  H2 
1590  PRINT 
1600  PRINT  *’ 

L  2  PH  2  TH1  TH2  PR  K" 

•410  ~RI NT  " 

>.  FT)  (DEG)  (DEG)  (DEG)  ( DB)  ( INT) " 

1  o20  Gl  =  HI  r  2 
.  >30  G2  =  H2  »  2 
1  ‘  40  G3  =  Gl  t  2 
*.)  G4  =  G2  t  2 
1650  G9  =  Gl  *  G2 
1  o  70  X  =  0 

1630  F0R  J  =  0  T0  J9 
159  0  K  =  0 

.700  M2  =  (LSCJ))  t  2 
i  /  10  M4  -  M2  t  2 
1  720  M 7  =  M2  *  Gl 
1730  M '  =  M2  *  G2 

1740  A4  =  K4  *  (  (M2+G1  +  G2)T2  -  4*G9  ) 

1750  A3  =  -2  *  K2  *  (  Q 1*  (M  4-  G9 )  +  G3  +  K2+G4  ♦  02*M7  +  Q3+M8  ) 
1 7 o L  A2  =  Q4+M4  +  G3  +  K4*G4  +  Q6+M7  +  97*M8  -  Q5*G9 

17  /0  A 1  =  -  2  *  (  Q  1*M  4  +  M7  ♦  K2*M8  ) 

1780  B4  =  4  *  A4 

1790  B3  =  3  *  A3 
1800  B2  =  2  *  A2 

18  10X2=  X  t  2 
\biO  X3  =  X2  *  X 
1830  X 4  =  X 3  *  X 


A-  3 


1840  F  =  A4*X4  +  A3*X3  +  A2*X2  +  A1*X  +  M4 

1850  G  =*  B4*X 3  +  B3*X2  +  B2*X  +  A1 

1860  E  =  F  /  G 

1870  X  =  X  -  E 

1880  K  =  K  +  1 

1390  IF  K  >  50  THEN  19  10 

1900  IF  ABSCE)  >  0.000001+X  THEN  18  10 

19  10  T 1  =  ATANC  SQRTC  X/C  1-X)  )  ) 

1920  S2  =  CK1  *  CT  NC  T 1 )  )  t  2 

1930  IF  S2  >=  1  THEN  S2  =  0.9999999999999 

19  40  T2  =  ATANC  SQRTC  S2/C1-S2)  )  ) 

1950  Cl  =  C0SCT1) 

1960  B  =  C0SCT2)  /  Cl 
1970  P7  =  2  *  R1  *  Cl 
1980  P2  =  1  +  B*Z  1 
1990  A  =  K  1  *  H2  /  B 
2000  P3  =  H  1  +  A 
2010  P4  =  HI  +  A  /  C  8*  2) 

2020  P9  =  ABSC  P8*P7/C  P2*  SQRTC  P3*P4)  )  > 

2030  IF  P9  <  IE- 10  THEN  P9  =  IE- 10 
2040  PSCI*J)  =  20  *  LOG  1 0C  P9  ) 

2050  T$C I  *  J)  =  C 18  0/PI )  *  T2 
2060  T 1  =  C  180/PI  )  *  T 1 
2070  S  =  " 

ZZZZZBBZZZZ. ZZXBBZZZZ. %ZZ8BZZ%%. ZZXBBZZZZ. ZBBZZZB/” 

2030  PRINT  IN  FORM  S:  LSCJ>*  FSCJ>*  Tl*  T$CI*J>*  P$CI*J>*  K 
209  0  NEXT  J 
2100  PRINT 
2110  NEXT  I 

2  120  OPEN  /ADATA/*  BINARY  OUTPUT*  1 
2130  OPEN  /CDATA/*  BINARY  OUTPUT*  3 

2140  WRITE  ON  1:  T$CI*J>  FOR  J  a  0  TO  J9  FOR  I  *  1  TO  10 
2150  WRITE  ON  3s  PSCI*J>  FOR  J  =  0  TO  J9  FOR  I  =  1  TO  10 
2160  CLOSE  1*  3 
2170  K9  =  0 

2180  OPEN  /UY DATA/*  BINARY  OUTPUT*  9 

1190  WRITE  ON  9:  PI*  Rl*  HI*  Dl*  U*  Y*  SI*  PO*  K9 

2200  CLOSE  9 
2210  PRINT 
2220  PRINT 

2230  PRINT  "THIS  COMPLETES  THE  OUTPUT  OF  DATA  FROM  PROGRAM  /USF/ 
1240  PRINT 

2250  PRINT  "OUTPUT  FROM  /USF/  IS  STORED  IN  /ADATA/  AND  /CDATA/." 

2260  PRINT 

2270  G2  =  TIME/60 

2  '30  INTEGER  G3 

2290  G3  =  G2  -  GS 

2300  PRINT  G3l  "  TOTPL  ELAPSED  SECONDS" 


PROGRAM  /AA/ 


1000  PRINT  "THIS  IS  PROGRAM  /AA/,  F0R  COMPUTING  THE  ANGLE  0F 
1010  PRINT  "  ARRIVAL  0F  UNDERWATER  SIGNAL  RELATIVE  T0  THE” 
1020  PRINT  *’  HORIZONTAL  AXIS  OF  A  SONAR  RECEIVER  AT  VARIOUS 
1030  PRINT  "  HEADINGS  AND  AT  VARIOUS  DEPTHS" 

10^0  PRINT 

1050  PRINT  "DATE  AND  TIME:  ":  DATE 

1060  PRINT 

1070  PRINT 

1080  VAR  =  ZERO 

1090  G 1  =  TIME/60 

1100  DIM  F$( -  60: 60) 

1110  REAL  T$(  l:  10,-60:  60),  A$C  1 :  3,  1 :  4,  -  60:  60> 

1120  INTEGER  PI,  Rl,  HI,  Dl,  U,  Y,  SS,  PO,  K9 

1130  IF  D 1  =  0  THEN  D 1  =  1  E-  6 

1140  OPEN  /UY DATA/,  BINARY  INPUT,  9 

1150  INPUT  FROM  9:  PI,  Rl,  HI,  Dl,  U,  Y,  S$,  PO,  K9 

1  160  CLOSE  9 

1170  II  =  0 

1 180  J9  =  Y  /  U 

1190  FOR  J  =  - J9  TO  J9 

1200  D2  =  J  *  U 

1^10  IF  Dl  =  0  AND  D2  <  0  THEN  FSCJ)  =  -90 

1220  IF  Dl  =  0  AND  D2  =  0  THEN  F$<J>  =  0 

1230  IF  Dl  =  0  AND  D2  >  0  THEN  F$(J)  =  90 

1240  IF  Dl  >  0  THEN  FJ(J)  =  (180/PI)  *  ATANC  D2/D1  ) 

1250  NEXT  J 

1260  OPEN  /ADATA/,  BINARY  INPUT,  1 

1270  INPUT  FROM  1:  T$(I,J)  FOR  J  =  0  TO  J9  FOR  I  =  1  TO  10 
1280  CLOSE  1 

1290  SO  =  "  *  DEPTH  OF  RECEIVER  (FEET)  =  •  XXXX//" 

1300  SI  =  "’HEADING  OF  RECEIVER  (DEGREES)  =’XXX//" 

1310  S2  =  " 1 0 ( XXX, XB) /" 

1320  S3  =  "  XXX. X B  /" 

1330  FOR  H2  =  200,  400,  800 

1340  11=  11+  1 

1350  PRINT  IN  FORM  SO:  H2 

1360  I  =  H2/U 

1370  PRINT 

1380  PRINT  " 

ANGLES,  THETA(A),  OF  SIGNAL  ARRIVAL  RELATIVE  TO  HORIZONTAL" 


1390  PRINT  " 

PERPENDICULAR  BISECTOR  0F  AXIS  0F  RECEIVER  C  DEGREES)  t  ** 

1400  FOR  K  =  1  TO  4 
1410  F2  =  30  *  CK-  1) 

1420  FOR  J1  =  1  TO  J9 

1430  J  =  -J 1 

1440  TSC I j  J)  =  T$C I  >  J  1 ) 

1450  NEXT  J1 

1460  FOR  J  =  - J9  TO  J9 

1470  A1  =  CPI/180)  *  C  F2  -  FSCJ)  ) 

1480  A2  =  CPI/180)  *  T$CI,J) 

149  0  B 1  =  C0SCA1)  *  SINCA2) 

1500  B2  =  SGNCB1) 

1510  83  =  B1  t  2 

1520  A$C  I  I*  K>  J)  =  B2  *  C  180/PI)  *  ATANC  SQRTC  B3/C  1-  B3)  )  ) 

1530  NEXT  J 
1540  PRINT 

1550  PRINT  IN  FORM  SI:  F2 

1560  PRINT  IN  FORM  S2:  A$CI1>K>J)  FOR  J  =  -J9  TO  -1 
1570  PRINT  IN  FORM  S3:  ASCI1,K>0) 

1580  PRINT  IN  FORM  S2:  A$(I1«K»J>  FOR  J  =  1  T0  J9 

159  0  PRINT 

1600  NEXT  K 

16  10  PRINT 

1620  NEXT  H 2 

1630  OPEN  /BDATA/,  BINARY  0UTPUT>  2 
1640  FOR  11=  1  TO  3 

1650  WRITE  ON  2:  ASCI1«K*J)  FOR  J  =  -J9  TO  J9  FOR  K  =  1  TO  4 
1660  NEXT  I  1 
1670  CLOSE  2 
i 68 0  PRINT 
1690  PRINT 

1700  PRINT  “THIS  COMPLETES  THE  OUTPUT  OF  DATA  FROM  PROGRAM  /AA/. 
1710  PRINT 

1720  PRINT  "OUTPUT  FROM  /AA/  IS  STORED  IN  /BDATA/." 

1730  PRINT 
1740  G2  =  TIME/60 
1750  INTEGER  G3 
1760  G3  =  G2  -  G 1 

1770  PRINT  G3:  "  TOTAL  ELAPSED  SECONDS" 


PROGRAM  /APL0T/ 


1000  PRINT  “THIS  IS  PROGRAM  /APLOT/,  FOR  PLOTTING  TIME  HISTORIES" 
1010  PRINT  "  OF  ANGLES  OF  ARRIVE.  AT  SONAR  RECEIVER" 

1020  PRINT 

1030  PRINT  "DATE  AND  TIME:  DATE 

1040  PRINT 

10S0  PRINT 

1060  VAR  =  ZERO 

1070  G 1  =  TIME/60 

1080  REAL  A$(  1:  3,  1 :  4,-60:  60) 

!J‘>)  INTEGER  PI,  Rl,  HI,  Dl,  U,  Y,  SI,  PO,  K9 
UCO  OPEN  /UY DATA/,  BINARY  INPUT,  9 
1  1  10  INPUT  FROM  9:  PI,  Rl,  HI,  Dl,  U,  Y,  SI,  PO,  K9 
1  120  CLOSE  9 
1 1 30  J9  =  Y  /  U 

1140  PRINT  "WHAT  IS  AIRCRAFT  SPEED,  Si,  IN  KNOTS": 

1  ISO  INPUT  SI 
1160  PRINT 

1170  TO  =  (3600/6076)  *  (U/Sl) 

1  18  0  T 1  =  20*  TO 

1190  OPEN  /BDATA/,  BINARY  INPUT,  2 
1200  OPEN  /APDATA/,  OUTPUT,  8 
1210  FOR  11=  1  TO  3 
1220  XO  =  500  +  (I  1-  D*  1000 
1230  YO  =  500 

1240  INPUT  FROM  2:  A$(I1,K,J)  FOR  J  =  -J9  TO  J9  FOR  K  =  1  TO  4 
1250  0  =  "WZZ%%3Z7.17.BZ/" 

1260  FOR  K  =  1  TO  4 
1270  Z  =  0 

1230  FOR  J  =  - J9  TO  J9 

1290  X  =  XO  +  J*T  1 

1300  Y  =  YO  +  2.  5*A$(  I  1,K,  J) 

1310  WRITE  ON  8  IN  FORM  0:  X,  Y,  Z 
1320  Z  =  1 
1330  NEXT  J 
1340  NEXT  K 

1350  NEXT  II 
1360  X  =  3500 
1370  Y  =  0 
138  0  Z  =  0 

1390  WRITE  ON  8  IN  FORM  0:  X,  Y,  7. 
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1400  CLOSE  8 
1410  CLOSE  2 
1420  PRINT 
1430  PRINT 

1440  PRINT  "THIS  COMPLETES  THE  OUTPUT  OF  DATA  FROM  PROGRAM  /APLOT/. 
1450  PRINT 

1460  PRINT  "OUTPUT  FROM  /APLOT/  IS  STORED  IN  /APDATA/> " 

1470  PRINT  "  READY  FOR  PUNCHING  OUT  ON  PAPER  TAPE.” 

1480  PRINT 
1490  G2  =  TIME/60 
1500  INTEGER  G3 
1510  G3  =  G2  -  G 1 

1520  PRINT  G3:  "  TOTAL  ELAPSED  SECONDS" 


PROGRAM  /ADDAMB/ 


1000  PRINT  "THIS  IS  PROGRAM  /ADDAMB/,  FOR  COMPUTING  ADDED  EFFECTS  OF' 
1010  PRINT  "  AMBIENT  NOISE  ON  TOTAL  UNDERWATER  SOUND  FIELD" 

1020  PRINT 

1030  PRINT  "DATE  AND  TIME:  ":  DATE 

1040  PRINT 

1050  PRINT 

1060  VAR  =  ZERO 

1070  G 1  =  TIME/60 

1080  REAL  P$(  It  10>0:60) 

1090  INTEGER  PI,  Rl,  HI,  Dl,  U,  Y,  SI,  PO,  K9 

1100  OPEN  /UY DATA/,  BINARY  INPUT,  9 

1110  INPUT  FROM  9:  PI,  Rl,  HI,  Dl,  U,  Y,  SI,  PO,  K9 

1  120  CLOSE  9 

1 130  J9  =  Y  /  U 

1140  PRINT  "WHAT  IS  AMBIENT  SOUND  PRESSURE  LEVEL,  PO,  IN  DB//UBAR": 
1150  INPUT  PO 
1160  PRINT 

1170  OPEN  /UYDATA/,  BINARY  OUTPUT,  9 

1180  WRITE  ON  9:  PI,  Rl,  HI,  Dl,  U,  Y,  SI,  PO,  K9 

1  190  CLOSE  9 

1200  P8  =  10»(P0/10) 

1210  OPEN  /CDATA/,  BINARY  INPUT,  3 

1220  INPUT  FROM  3:  PS(I,J>  FOR  J  =*  0  TO  J9  FOR  I  =  1  TO  10 

1230  CLOSE  3 

1240  FOR  I  =  1  TO  10 

1250  FOR  J  =  0  TO  J9 

1260  P 7  =  10t(P$(I, J)/10) 

1270  P$C1,J)  =  1 O+LOG 1 0( P7+P8 ) 

1280  NEXT  J 
1290  NEXT  1 

1300  OPEN  /FDATA/,  BINARY  OUTPUT,  6 

1310  WRITE  ON  6:  PS<I,J)  FOR  J  =  0  TO  J9  FOR  I  =  1  TO  10 
1320  CLOSE  6 
1330  PRINT 
1340  PRINT 

1350  PRINT  "THIS  COMPLETES  THE  OUTPUT  OF  DATA  FROM  PROGRAM  /ADDAMB/. 
1360  PRINT 

1370  PRINT  "OUTPUT  FROM  /ADDAMB/  IS  STORED  IN  /FDATA/." 

1380  PRINT 
1390  G2  =  TIME/60 
1400  INTEGER  G3 
1410  G3  =  G2  -  G 1 

1420  PRINT  G3:  "  TOTAL  ELAPSED  SECONDS" 


PROGRAM  /LPL0T/ 


1000  PRINT  "THIS  IS  PROGRAM  /LPL0T//  P0R  PLOTTING  TIME  HISTORIES" 
1010  PRINT  "  OF  USF  SPLS  AT  SONAR  RECEIVERS  AT  VARIOUS  DEPTHS" 
1020  PRINT 

1030  PRINT  "DATE  AND  TIME:  "J  DATE 

1040  PRINT 

1050  PRINT 

1060  VAR  =  ZERO 

1070  G 1  =  TIME/60 

1080  REAL  P$< It  1 0/  0: 60) 

1090  INTEGER  PI/  Rl/  HI*  Dl*  U*  Yz  SI/  PO/  K9 
1100  OPEN  /UYDATA/z  BINARY  INPUT/  9 
1110  INPUT  FROM  9:  PI/  Rl*  HI/  Dl/  U/  Y/  SI/  PO/  K9 
1120  CLOSE  9 

1130  PRINT  "WHAT  IS  AIRCRAFT  SPEED/  SI/  IN  KNOTS": 

1140  INPUT  SI 
1150  PRINT 
1160  J9  =  Y  /  U 

1170  TO  =  (3600/6076)  *  (U/Sl) 

1180  T 1  =  20* TO 

1190  PRINT  "AMBIENT  SOUND  PRESSURE  LEVEL  IS  ":  PO:  "  DB//UBAR" 
1200  PRINT 

1210  OPEN  /FDATA/z  BINARY  INPUT/  6 

1220  INPUT  FROM  6:  PSCI/J)  FOR  J  =  0  TO  J9  FOR  I  =  1  TO  10 
1230  CLOSE  6 

1240  OPEN  /LPDATA/,  OUTPUT/  8 
1250  0  =  "ZZZZZXBZZZZBZ/” 

1260  XO  =  500 

1270  YO  =  500 

1280  FOR  I  =  1/  2/  4/  8 

1290  Z  =  0 

1300  FOR  J  =  - J9  TO  J9 

1310  X  =  XO  +  J*T1 

1320  Y  =  YO  +  5*PS(  I /  ABS(  J)  ) 

1330  WRITE  ON  8  IN  FORM  0:  X.*  Y/  Z 

1340  Z  =  1 

1350  NEXT  J 

1360  NEXT  I 

1370  Y  =  YO  +  5*P0 

1380  FOR  J  =  -J9  TO  J9  BY  6 

1390  X  =  XO  +  J* T 1 
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1400  Z  =  0 

1410  WRITE  0N  8  IN  FORM  0:  X,  Y,  Z 
1420  X  =  X  +  4*T1 
1430  Z  =  1 

1440  WRITE  0N  8  IN  FORM  0i  X,  Y,  Z 

1450  NEXT  J 

1460  X  =  1500 

1470  Y  =  0 

1480  Z  =  0 

'490  WRITE  0N  8  IN  FORM  0:  X>  Y>  Z 

i  >00  CLOSE  8 
1510  PRINT 
1520  PRINT 

1530  PRINT  "THIS  COMPLETES  THE  OUTPUT  OF  DATA  FROM  PROGRAM  /LPL0T/." 
1540  PRINT 

1550  PRINT  "OUTPUT  FROM  /L  PL  0  T  /  IS  STORED  ON  /LPDATA/,  " 

1560  PRINT  "  READY  FOR  PUNCHING  OUT  ON  PAPER  TAPE." 

1570  PRINT 
1580  G2  =  TIME/60 
1590  INTEGER  G3 
1600  G3  =  G2  -  G 1 

1610  PRINT  G3J  "  TOTAL  ELAPSED  SECONDS" 
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PROGRAM  /CON/ 

i 


1000  PRINT  “THIS  IS  PROGRAM  /CON/,  FOR  COMPUTING  EQUI PRESSURE  CONTOURS” 
1010  PRINT  "  ASSOCIATED  WITH  UNDERWATER  SOUND  FIELD  INCLUDING  AMBIENT** 
1020  PRINT 

1030  PRINT  "DATE  AND  TIME:  **:  DATE 

1040  PRINT 

1050  PRINT 

1060  VAR  =  ZERO 

1070  G 1  =  TIME/60 

1080  INTEGER  PI,  Rl,  HI,  Dl,  U,  Y,  S$,  PO,  K9,  Y2 
1090  INTEGER  L$C  1:  50) 

1100  REAL  PSC l: 10, 0: 60),  ESC  1:50) 

1110  REAL  YSC  1:200),  Z$C!?200) 

1120  OPEN  /UYDATA/,  BINARY  INPUT,  9 

1130  INPUT  FROM  9:  PI,  Rl,  HI,  Dl,  U,  Y,  S$,  PO,  K9 
1  140  CLOSE  9 
1150  J9  =  Y  /  U 

1160  OPEN  /FDATA/,  BINARY  INPUT,  6 

1170  INPUT  FROM  6:  PSCI,J)  FOR  J  =  0  TO  J9  FOR  I  =  1  TO  10 
1180  CLOSE  6 

1190  PRINT  "WHAT  IS  MAXIMUM  DISTANCE,  Y2,  TO  WHICH  UNDERWATER** 

1200  PRINT  "  EQUIPRESSURE  CONTOURS  ARE  TO  BE  COMPUTED,  IN  FEET": 

1210  INPUT  Y 2 
1220  PRINT 
1230  J9  =  Y2/U 

1240  PRINT  "THE  MAXIMUM  UNDERWATER  SPL  IS  **:  I NTC  PSC  1, 0) )  1  *’  DB//UBAR” 
1250  PRINT 

1260  PRINT  "THE  MINIMUM  UNDERWATER  SPL  IS  ”:  I  NTC  PSC  1,  J9  >  )  :  *’  DB//U3AR" 
1270  PRINT 

1280  PRINT  "WHAT  IS  DESIRED  CONTOUR  INTERVAL,  C2,  IN  THE" 

1290  PRINT  **  SPECIFIED  UNDERWATER  REGION,  IN  DB": 

1300  INPUT  C2 
1310  PRINT 

1320  OPEN  /EDATA/,  BINARY  OUTPUT,  5 
1330  DEF  FNACP)  =  10  t  CP/20) 

1340  K  1  =  I  NTC  PSC  1,  0)  /  C2  ) 

1350  K  2  =  I  NTC  P$C1,J9)  /  C2  ) 

1360  K9  =  K1  -  K 2 

1 370  PRINT  "THERE  ARE  **:  K9:  **  EQUI  PRESSURE  CONTOURS  IN  THE” 

1380  PRINT  "  SPECIFIED  UNDERWATER  REGION" 

139  0  PRINT 
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1400  FOR  K  =  1  TO  K9 
1410’tIStKT"::  *0*  *  ~ 

1420  ESCK)  =  CK1+1-K)  *  C2 

1430  NEXT  K 

1440  FOR  K  =  1  TO  K9 

1450  FOR  J  =  0  TO  J9 

1460  FOR  I  =  1  TO  9 

1470  F2  =  P$CI,J)  -  P$( I + 1 > J) 

1480  F  =  C  PSCI/J)  -  ES(K)  )  /  F2 
1490  IF  F  <  0  OR  F  >=  1  THEN  1570 
1500  LSCK)  =  LSCK)  +  1 
1510  L  =  LSCK) 

1520  FI  =  FNAC  PSCI>J)  )  -  FNAC  ESCK)  ) 

1530  F2  =  FNAC  PSC1>J)  )  -  FNAC  P$CI+1*J)  ) 

1540  F  =  F1/F2 

1550  YSCL)  =  J  *  U 

1560  ZSCL)  =  C  I  +  F  )  *  U 

1570  NEXT  I 

1580  NEXT  J 

1590  FOR  I  =  1  TO  10 

1600  FOR  J  =  0  TO  C J9-  1) 

1610  F2  =  PSCI, J)  -  PSCI,J+1) 

1620  F  =  C  PSCI/J)  -  ESCK)  )  /  F2 
1630  IF  F  <  0  OR  F  >=  1  THEN  1710 
1640  LSCK)  =  LSCK)  +  1 
1650  L  =  LSCK) 

1660  FI  =  FNAC  PSCI,J)  )  -  FNAC  ESCK)  ) 

1670  F2  =  FNAC  PSCI,J)  )  -  FNAC  P$CI,J+1)  ) 

1680  F  =  F1/F2 

1690  YSCL)  =  C  J+F  )  *  U 

1700  ZSCL)  =  I  *  U 

1710  NEXT  J 

1720  NEXT  I 

1730  FOR  I  =  1  TO  9 

1740  FOR  J  =  0  TO  CJ9-  1) 

1750  F2  =  PSCI,J)  -  P$CI+1,J+1) 

1760  F  =  C  PSCI/J)  -  ESCK)  )  /  F2 
1770  IF  F  <  0  OR  F  >=  1  THEN  1850 
1780  LSCK)  =  LSCK)  +  1 
1790  L  =  LSCK) 

1800  FI  =  FNAC  PSC  I»  J)  )  -  FNAC  ESCK)  ) 

1810  F2  =  FNAC  P$CI,J)  )  -  FNAC  P$CI  +  1,J+1)  ) 

1820  F  =  F1/F2 

1830  YSCL)  =  C  J+F  )  *  U 

1840  ZSCL)  =  C  I  +  F  >  *  U 

1850  NEXT  J 

1860  NEXT  I 
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1870  FOR  I  =  1  TO  9 

188 0  FOR  J  *  1  TO  J9 

1890  F2  =  P$CI#J)  -  P$CI  +  1#J-1) 

1900  F  =  C  P$CI#J)  -  E$<K)  )  /  F2 
19  10  IF  v  <  0  OR  F  >=  1  THEN  1990 
1920  LSCK)  =  L$<K)  +  1 
1930  L  =  LSCK) 

1940  FI  *  FNAC  P$CI#J)  )  -  FNAC  ESCK)  ) 

1950  F2  =  FNA<  P$CI#J)  )  -  FNAC  P$CI+1#J-1)  ) 

1960  F  =  F1/F2 

1970  Y $CL)  =  C  J-F  )  *  U 

1980  Z$(L)  =  <  I  +  F  )  *  U 

1990  NEXT  J 

2000  NEXT  I 

2010  SI  »  "ZXXXB" 

2020  PRINT  IN  FORM  SI:  ESCK) 

2030  PRINT  "D8  CONTOUR  POINTS’* 

2040  FOR  LI  =  1  TO  LSCK) 

2050  Z9  =  -  1 

2060  FOR  L  =  1  TO  LSCK) 

2070  IF  ZSCL)  <  Z9  THEN  2100 
2080  Z9  =  ZSCL) 

209  0  L9  =  L 

2100  NEXT  L 

2110  S2  =  "XXXXX.XX BBXXXXX.XZ/" 

2120  PRINT  IN  FORM  S2:  YSCL9)#  ZSCL9) 

2130  WRITE  ON  5:  YSCL9)#  ZSCL9) 

2140  ZSCL9)  »  -1 
2150  NEXT  LI 
2160  PRINT 
2170  NEXT  K 

2180  OPEN  /DDATA/#  BINARY  OUTPUT#  4 
2190  WRITE  ON  4:  LSCK)  FOR  K  =  1  TO  K9 
2200  CLOSE  4 
2210  CLOSE  5 

2220  OPEN  /UYDATA/,  BINARY  OUTPUT#  9 

2230  WRITE  ON  9:  PI#  Rl#  HI#  Dl#  U#  Y#  S$#  PO#  K9#  Y2 
2240  CLOSE  9 
2250  PRINT 
2260  PRINT 

2270  PRINT  "THIS  COMPLETES  THE  OUTPUT  OF  DATA  FROM  PROGRAM  /CON/." 
2280  PRINT 

2290  PRINT  "OUTPUT  FROM  /CON/  IS  STORED  IN  /DDATA/  AND  /EDATA/." 

2300  PRINT 

2310  G2  *  TIME/60 

2320  INTEGER  G3 

2330  G3  =  G2  -  G1 

2340  PRINT  G3t  "  TOTAL  ELAPSED  SECONDS" 
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PROGRAM  /CPL0T/ 


1000  PRINT  "THIS  IS  PROGRAM  /CPL0T/,  FOR  PLOTTING  EQUI  PRESSURE** 
1010  PRINT  "  CONTOURS  FOR  USF  SPL  S  INCLUDING  AMBI  ENT  EFFECTS” 
1020  PRINT 

1030  PRINT  "DATE  AND  TIME:  DATE 

1 040  PRINT 

1050  PRINT 

1060  VAR  =  ZERO 

1070  G 1  =  TIME/60 

1080  INTEGER  L$<  1:  50) 

1090  REAL  YSC  1:200),  Z$<  1:200) 

1100  INTEGER  PI,  Rl,  HI,  Dl,  U,  Y,  SI,  PO,  K9,  Y2 
1  1  10  INTEGER  E 

1120  OPEN  /UYDATA/,  BINARY  INPUT,  9 

1130  INPUT  FROM  9:  PI,  Rl,  HI,  Dl,  U,  Y,  SI,  PO,  K9,  Y2 
1140  CLOSE  9 
1150  J9  =  Y2  /  U 
1160  C  =  50/U 

1170  OPEN  /DDATA/,  BINARY  INPUT,  4 
1180  INPUT  FROM  4:  L$(K)  FOR  K  =  1  TO  K9 
1190  CLOSE  4 

1200  OPEN  /EDATA/,  BINARY  INPUT,  5 
1210  OPEN  /CPDATA/,  OUTPUT,  8 
1220  0  =  "ZZZZZZBZZZZBZ/" 

1230  XI  =  0 
1240  Y 1  =  0 
1250  XO  =  0 
1260  YO  =  600 
1270  X  =  XO 
1280  Y  =  YO 
1290  Z  =  0 


1300 

WRITE  ON 

8  IN 

FORM 

0: 

V 

A/ 

Y, 

Z 

1310 

X  =  XO  + 

C*Y  2 

1320 

Z  =  1 

1330 

WRITE  ON 

8  IN 

FORM 

0: 

X, 

Y, 

Z 

1340  FOR  K  =  1  TO  K9 
1350  Z  =  0 

1360  FOR  L  =  1  TO  L$CK) 

1370  INPUT  FROM  5:  YS(L),  ZSCL) 
1380  NEXT  L 
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1390  FOR  L  =  1  T0  LSCK) 

1400  X  =  XO  +  C*YSCL) 

1410  Y  =  Y0  -  C*Z$<L> 

1420  IF  ABSCX-X 1)  <  50*C  AND  ABS(Y-Yl)  <  20*C  THEN  1470 
1430  XI  =  X 
1440  Y 1  =  Y 

1450  WRITE  0N  8  IN  F0RM  0:  X,  Y,  Z 
1460  Z  =  1 
1470  NEXT  L 
1430  NEXT  K 

1490  IF  H 1  >=  4* U  THEN  P9  =  PI  -  2O+L0GC  CH1-4*U)/R1  ) 

1500  IF  H1<4*U  THEN  P 9  =  PI  -  30 

lc10  P8  =  PI  -  2O*L0G1O<  (  SQRTCHlt2+(J9*U)  t2))/Rl  ) 

1520  PRINT  **THE  MAXIMUM  SPL  IN  THE  ATMOSPHERIC  REGION” 

1530  PRINT  ”  0F  INTEREST  IS  INT(P9>:  ”  DB//UBAR” 

1540  PRINT 

1550  PRINT  ”THE  MINIMUM  SPL  IN  THE  ..10SPHERI C  REGION” 

1560  PRINT  ”  OF  INTEREST  IS  ”:  INT(P8):  ”  DB//UBAR" 

1570  PRINT 

1580  PRINT  "WHAT  IS  DESIRED  C0NT0UR  INTERVAL,  Cl,  IN  AIR  JUST  ABOVE” 
159  0  PRINT  ”  SURFACE  OF  SPECIFIED  UNDERWATER  REGION,  IN  DB”: 

1600  INPUT  Cl 
16  10  PRINT 

1620  K 1  =  I NT(  P9/C1  ) 

1630  K2  =  I NT(  P8/C1  ) 

1640  K8  =  K1  -  K2 
.650  FOR  K  =  1  T0  K8 
'660  E  =  (K 1+  1-K)  *  Cl 

1670  PRINT  "CONTOUR  K:  ”  IN  AIR  IS  ”:  E:  ”  DB//UBAR” 

1680  R  =  C  *  R1  *  10U  <Pl-E)/20  ) 

1690  Z  =  0 

1700  FOR  L  =  0  TO  Y2  BY  10 

17 10  X  =  XO  +  C*L 

1720  IF  X  >  R  THEN  1790 

1730  Y  =  YO  +  C*H  1  -  S9RT<  Rt  2-X»2) 

1740  IF  Y  <=  YO  THEN  1730 

1750  WRITE  ON  8  IN  FORM  0:  X,  Y,  Z 

1760  IF  Y  >  YO  +  200  THEN  1790 

1770  Z  =  1 

1780  NEXT  L 

1790  NEXT  K 

1800  X  =  XO  +  C*Y 2  +  500 
18  10  Y  =  0 
1820  Z  =  0 

1530  WRITE  ON  8  IN  FORM  0:  X,  Y,  Z 

1340  CLOSE  8 

1350  CLOSE  5 

18  60  PRINT 

1370  PRINT 
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1880  PRINT  “THIS  COMPLETES  THE  OUTPUT  0F  DATA  FROM  PROGRAM  /C PLOT/. 
1890  PRINT 

1900  PRINT  "OUTPUT  FROM  /CPLOT/  IS  STORED  IN  /  CP  DATA/, " 

19  10  PRINT  "  READY  FOR  PUNCHING  OUT  ON  PAPER  TAPE." 

1920  PRINT 
1930  G2  =  TIME/60 
1940  INTEGER  G3 
19  50  G3  =  G2  -  G 1 

1960  PRINT  G3:  "  TOTAL  ELAPSED  SECONDS" 
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